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Abstract 

We study the possibility of identifying the CP violating phases in the PMNS mixing matrix in 
the lepton sector and also that in the CKM mixing matrix in the quark sector with the phase 
responsible for the spontaneous CP violation in the Higgs potential, and some implications. Since 
the phase in the CKM mixing matrix is determined by experimental data, the phase in the lepton 
sector is therefore also fixed. The mass matrix for neutrinos is constrained leading to constraints 
on the Jarlskog CP violating parameter J, and the effective mass < mpp > for neutrinoless double 
beta decay. The Yukawa couplings are also constrained. Different ways of identifying the phases 
have different predictions for \x — > eee and r — > hhh- Future experimental data can be used to 
distinguish different models. 
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I. INTRODUCTION 



The origin of CP violation is one of the most outstanding problems of modern particle 
physics. In the quark sector, the leading effect of CP violation comes from the phase 
(Jkm, which sometimes is referred as the Dirac phase, in the Cabbibo-Kobayashi-Maskawa 
(CKM) model jl, 2|. CP violation may also exist in the lepton sector. One of the ways to 
describe low energy CP violation in this sector is to have CP violating Dirac phase 5pmns m 
the Pontcove-Mkai-Nagawa-Sakata (PMNS) mixing matrix {3, 4]. If neutrinos are Majorana 
particles, besides the Dirac phase 5pmns ; there may also exist CP violating Majorana phases. 
It is important to understand the origin of CP violation. An interesting proposal due to 
T.-D. Lee is that CP is spontaneously violated 0]. The popular Weinberg model js] of 

8[ with data and has been decisively ruled 
. In a previous work, we have studied 



spontaneous CP violation model has problems 
out by CP violating measurement in B decays 



the possibility of restoring the idea that CP is broken spontaneously and identifying the 
phase 5km, up to a sign, as the phase S spon that causes spontaneous CP violation in the 
Higgs potential [lOj . In this work we generalize this approach to the lepton sector to study 
how CP violating phase <5pmns ma Y be identified with the phase S spon , and some of the 
consequences. 

A model to realize the above idea for quark sector has been studied in Ref. 10] . It 
involves multi-Higgs fields. In that work, in order to avoid possible strong CP problem a 



Peccei-Quinn symmetry is imposed to the model [8(|. Implementing spontaneous CP 



violation with PQ symmetry then requires more than two Higgs doublets [ijj]. For our 
purpose we find that in order to have spontaneous CP violation with PQ symmetry at least 



three Higgs doublets 4> 
singlet S = e i6s S = e { 



V2 



[v, + Ri + iAA, hi 



and one complex Higgs 



s + R s + iA s )/\/2 are needed. The Higgs singlet with a large 
vacuum expectation value (VEV) renders the axion from PQ symmetry breaking to be 
invisible {13I . Q] , thus satisfying experimental constraints on axion couplings to fermions. 



Two possible ways of assigning PQ charges to the quarks and Higgs bosons were consid- 



ered 



10| 



Model a) Q L : , U R : -1 , D R : -1 , lj2 : +1 
Model b) Q L : , U R : +1 , D R : +1 , <f> 1)2 : +1 



■i; 
-1 . 



(i) 



2 



In both cases, S has PQ charge +2. 

With these PQ charge assignments, it is possible to have spontaneous CP violation with 
a non-zero phase of VEVs, <5 spon = 9 2 — Q\. The other phases #3 and 6 S are related to <5 spon 
by minimization conditions. 

The resulting quark mass terms in the Lagrangian, after removing un-physical phases 
and the diagonalized basis of down quark mass matrix M d and up quark mass matrix M u 



for models a) and b), respectively, we have 10| 



Model a) : C m = -Ur 



M ul + M u2 e i5 ^ U R - D L M d D R + h.c. 



Model b) : C m = -U L M U U R - D L M dl + M d2 e 



20 spoil 



D R + h.c. , 



(2) 



where M ui = —Y ui Vi/\/2 and M di = —T di Vi/\/2 are all real since spontaneous CP violation 
is imposed. 

The above mass matrices can be transformed into diagonalized form M Ujd by bi-unitarity 
transformation. For Model a), M u = V CKM M u V R j , and for Model b), M d = V^ KM M d V R j . A 
direct identification of the phase <5 spon with the phase 5km m the CKM matrix is not possible 
in general at this level. There are, however, classes of mass matrices which allow such a 
connection. A simple example is provided by setting V R d to be the unit matrix. With this 



condition, one has 



io| 



Model a) : V^ KM = (M ul + e iSspou M u2 )M~ l 

Model b) : V^ckm = (M dl + e~ iSapon M d2 ) Mj 1 (3) 

Expressing the CKM matrix in the above form is very suggestive. If Vqkm (or ^ckm) can 
always be written as a sum of two terms with a relative phase, then the phase in the CKM 
matrix can be identified with the phase 5 sp0 n- 

Since multi-Higgs bosons are involved, there are in general flavor changing neutral current 
(FCNC) interactions mediated by neutral Higgs bosons at the tree level. However, unlike 
in general multi-Higgs models, the FCNC Yukawa couplings are fixed in terms of the quark 
masses, CKM mixing angles, and vacuum expectation values Vi of Higgs bosons making 
phenomenological studies much easier and having different predictions than general multi- 
Higgs models [k], [ijj] . In the following we will study how the same idea can be applied to 
lepton sector and related implications. 



II. CP VIOLATING PHASE IN THE LEPTON SECTOR 



There are some complications for a straightforward generalization of the idea discussed 
in the previous section to lepton sector because neutrinos may be Majorana type. Also in 
order to have non-zero neutrino masses in the model, some extensions are needed. We will 
introduce right-handed neutrinos to achieve this through seesaw mechanism. 

The charged current mixing matrix Vpmns i n the lepton sector is given by, 

£ = -X h-fVpuNsVLW- + h.c. (4) 

with Vpmns = VlYi}- Here V[ and V£ are unitary matrices transforming the charged lepton 
mass matrix Mi and the neutrino mass matrix M v to their diagonalized form M\ and M u , 
Mi = vSMiVjt and M v = V?M V V£*. 

We consider two possible ways given in the discussion of Ref. [h]] to couple the leptons 
to Higgs bosons parallel to the quark sector with PQ charges assigned to leptons as 

Model a): L L (0) , l R (-l) , N R (-1), 

Model b): L L (0), l R (+l) , N R (+1) . (5) 

Combining with PQ charges for Higgs in Eq. the Yukawa interaction of above two 
models is then given by [lo| 



Model a) : C = L L {Y X H X + Y 2 H 2 e i5 °»™)N R + L L Y 3 H 3 l R + (N R yY s SN R + h.c, 
Model b) : C = L^Y 3 H 3 N R + T L (J X H X + Y 2 H 2 e- i5 °*°»)l R + (NrYY^Nr + h.c. , (6) 

where Y x , Y 2 , Y 3 , and Y s are 3x3 real coupling matrices. The interaction shown above 



gives the mass terms of lepton sector as 10] 



C m = -l L Mil R - T7ZM D N R - - {N R YM R N R , (7) 



with Mo, Mi, and M R defined below 



Model a): Mi = -^=Y 3 v 3 , M D = -^=(Y lVl + Y 2 v 2 e t6 —) , M R = —y/2Y s v s , 
Model b): Mi = -^=(Y lVl + Y 2 v 2 e- 15 ^), M D = ~^=Y 3 v 3 , M R = -V2Y s v s . 
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To get neutrino masses, we expressed the neutrino mass terms in the usual seesaw mass 

n 

matrix with the same notation in Ref. [16[ shown as 



1 



h.c. 



V 



(9) 



where M seesaw is defined as 6 x 6 mixing matrix, with 



Mr 



D 



Mr 



(10) 



Strictly speaking V^} is only approximately unitary. The relation between Mr, and di- 



agona 



by |1Q(| 



ized light neutrino mass matrix M„, also with heavy one M^ = Mr 16l. 117]. is given 



M D M^M T D = -V?M V VZ* = (iV?)M y (iVr)- 



The general solutions for Mp are given by 18] 



M D = (iV?)Ml/ 2 OMH 2 . 



(11) 



(12) 



where O is a matrix satisfying 00 T = I. Choosing O = I is similar to choosing Vr = I for 
Model b) in the quark sector. 

We now discuss several simple ways to identify phases in the mass matrix with the 
spontaneous CP violating phase 5 spon . 



A. Phase in Model a) 

For Model a), we work with the basis where Mi is already diagonalized. If now we redefine 
II to cause V[ as i, then we have Vpmns = iVj} and one can get rid of the '-' sign in the 
Eq. flTJD and set 

M D = M m + M D2 e^°» = V PMm M l J 2 OM]l 2 . (13) 

where Mm = —Y(Vi/ y2 are real matrices. Note that for Ir we also have V R = i at the same 
time. 
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In general one can express Mm,2 as follows 

M m = ReiVpMNsM^OMll 2 ) - cot 5 spon lm(V PM m 2 OM^ 2 ) , 

M m = — | \m{V ?M ^Ml' 2 OM l J 2 ) . (14) 

sin S p On 

There are some complications compared with realization of the same in the quark sector 
because the Majorana feature of neutrinos. The Vpmns matrix can be separated into two 
parts Vpmns — Vpmns V^>. Here Vpmns is a unitary matrix with only Dirac type of phase. The 
diagonal phase matrix V p can be chose by convention to be V p = diag(l, e mi , e l<72 ). These 
phases can only be measured in AL = 2 processes such as neutrinoless double (3 decay. 
There is also the possibility that phases appear in the matrix O. 

Identifying S spon with the phase in Vpmns is the emphasis of the text. The phase, in 
general, can appear as the Dirac phase in Vpmns- If the Vpmns does not have CP violation, 
for example it is a tri-bimaximal form, then one has to identify S spon with the phases in V p 
or in O. We will concentrate on the possibility of identifying <5 spon with the Dirac phase in 

VpMNS- 

In this case, the spontaneous CP violating phase <5 S pon is identified with the Dirac phase 
in Vpmns which is the only phase in the model. For the mixing matrix, we can work with the 
Particle Data Group (PDG) parametrization of the mixing matrix [9], but with modification 
to write it in a form with just one phase appearing with the same sign for simple identification 



with 5 spon 1 1 



yPDG 
^PMNS 



( ci 2 c 13 e iS ™ s 12 c 13 e i5 ^ s l3 N ^ 

-si 2 c 23 - c l2 s 23 si 3 e thi C12C23 - s 12 S2 3 Si 3 e l5w s 23 c 13 
\ si 2 s 23 - c 12 c 23 s 13 e lSl3 -c 12 s 23 - s 12 c 23 si 3 e tSvi c 23 c 13 ) 



(15) 



where Sy = sin^- and c^- = cos 6^-. 

The matrices with above parametrization are related to the usual PDG parametrization 
by multiplying a diagonal phase matrix Vs = diag(e~ l<513 , 1, 1) on left-handed side. We make 
the Dirac phase 6i 3 in Vpmns to be identical to 5 spon , and also because 5 sp on is identical with 
Dirac phase in CKM matrix in our model, we take the global fitting results from PDG 9] 
to get 5i 3 = (68.9 ± 2.7)° for CKM matrix with modified PDG parametrization. We have 
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the mass matrices Mdi, 2 given by 

/ 



M 



PDG 
Dl 



\ 



M}/ 2 OMl /2 



L N 



(16) 



M 



PDG 
D2 



M}J 2 0m1 12 



(17) 



s 13 
-S12C23 C12C23 S23C13 
\ S12S23 -C12S23 C23C13 J 
1 C12C13 S12C13 » 
-C12S23S13 -S12S23S13 

V -C12C23S13 -S12C23S13 j 

The non-zero Dirac 5x3 implies none of the mixing angles in the Vpmns can be zero. 

Note that the solutions discussed above are not unique. This implies that in order to 
completely fix a underlying model, more physical requirements are needed. To see the 
non-uniqueness, we take another parametrization for the mixing matrix, the original KM 
matrix [2|, for illustration. There the mixing matrix is parameterized as 



KM 
PMNS 



-S1C3 



-S1S3 



sic 2 C1C2C3 - s 2 s 3 e lS C1C2S3 + s 2 c 3 e lS 
\ sis 2 C1S2C3 + c 2 s 3 e lS C1S2S3 - c 2 c 3 e lS J 



(18) 



with Si = sin 9i and q = cos^j. In this case one identifies 5 = 5 spon with 6 = (90.3 ± 2.7)° 



^ ci -S1C3 -S1S3 ^ 



obtained from PDG global fitting results [9|, and the mass matrices M D1 2 can be written as 

Ml^OMH 2 , (19) 

3/ 



M 



KM 
Dl 



M KM 
1VI D2 



m}J 2 om]1 2 



(20) 



sic 2 C1C2C3 C1C2S3 

\ s l s 2 CiS 2 C 3 C\S 2 S 3 J 

/ 

-s 2 s 3 s 2 c 3 
y c 2 s 3 -c 2 c 3 J 

We emphases that the above two ways of parameterizing the mixing matrices are two 
different models. We will refer them as Model a(PDG)) and Model a(KM)). The differences 

n 

between these two models will show up in the CP violating Jarlskog parameter [19J J, in 
neutrinoless double beta decay and also in the Yukawa couplings since Mm for these two 
models are different. 

Since we have identified the Dirac phase with the spontaneous CP violating phase which 
are also related to the quark mixing CP violating phase, the Dirac phase in Vpmns is known. 



For Model a(PDG)), using the current range of mixing parameters determined from various 
neutrino oscillation data |20| 

0x2 = 34.5 ± 1.0 (tli)° , 9 23 = 42.8^; 9 7 (^ 3 7 )° , 9 13 = 5.1j£°(< 12.0)°, (21) 

we obtain the central and la allowed range for J(PDG) 

J(PDG) = 0.019 ±0.012. (22) 

For Model a(KM)), we can obtain the central and la ranges for the mixing angles from 
I Kill IK2I, and |V e 3| by using the values in Eq. f[2"Tj) . The results are 

sin 9 1 = 0.57 ±0.02; sin 9 3 = 0.16 ± 0.10. (23) 

Since with this model one also identifies the phase in Eq. ( fl8|) with Dirac phase in the 
original KM parametrization of quark sector, that is, 5 = (90.3 ± 2.7)° as shown before, 
Using this value with Eq. (1231) we can derived the corresponding values for s 2 from {V^l or 
\V T 3\, which is given by 

sin 9 2 = 0.68 ±0.04. (24) 

The value for Jarlskog parameter in this case J(KM) is different than that for J(PDG). 
Using the above derived values for mixing angles and phase, we have 

J (KM) = 0.021 ±0.013. (25) 

In principle measurement of J can be used to distinguish different models, but the difference 
is small to it practically difficult. 

The above two models also have difference in predicting the effective mass < 171/3/3 >= 
I J2i m iVei\ f° r neutrinoless double beta decay which we discuss in the following. 

At present, the absolute mass scale for neutrino is not known although stringent constraint 



exists from cosmological consideration, which gives [21] J2i m i < 0.28eV. Combining with 
data from neutrino oscillation 



Am 2 12 = (7.59 ±0.20) x KT 5 (eV) 2 [22j; (26) 
Am 2 23 = (2.43 ±0.13) x 10" 3 (eV) 2 [23j, (27) 



one can constrain < mpp >. 
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0.00 [ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ' — 

0.00 0.02 0.04 0.06 0.08 



FIG. 1: < rripp > vs. mi for model a), with normal hierarchy, respectively. 

We will take the central values of the mixing angles and phases, neutrino mass squared 
differences, and also take neutrino mass m\ as a free parameter satisfying the cosmological 
constraint for illustration. The results are show in Figure [TJ We have only plotted the case 
with normal neutrino mass hierarchy. The results with inverted hierarchy are almost the 
same and have not been plotted. As can be seen from the figure that the two different models 
lead to similar results. There the < mpp > are in principle different but the differences are 
small. 

Another quantity related to neutrino mass is < m„ e >= (X)jWi 2 I^eil 2 ) 1 ^ 2 which can be 
measured in tritium decay. This quantity is however, the same for the two models we 
considered above. With the central values for mixing angles, it is given by 0.088 for normal 
hierarchy with vii\ = 0.088 close to its upper bound, for example. For inverted hierarchy 
with the same m 1; we get almost the same results as that we have in normal hierarchy. This 
quantity is not useful in distinguishing different models considered above. 

B. Phase in Model b) 

For Model b), information about CP violation is encoded in Mi. Mo is a real matrix 
from Yukawa coupling Y" 3 . If one recalls Eq. ( Tl2l) . it is straightforward to deduce that VI 
must be complex. For simplicity we will work in a basis with V£ = i and we have a simpler 
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relation 



M D = Ml' 2 OM]l\ 



(28) 



where O is a complex orthogonal matrix as mentioned before, we then follow the same 
philosophy to have set a possible non-unit unitary matrix Vr to i in order to getting the 
useful relation which is given as 



Mi = M a + M l2 e i5 ^ on = Vj MNS M t 



(29) 



In this basis Vpmns — ~^[, and we can apply the same idea in model (a) that Dirac phase 
Vpmns is identical to 5 spon . We also work with two models with the PDG and original KM 
parameterizations, Model b(PDG)) and Model b(KM)). For Model b(PDG)), we have 



Mi 



PDG 
11 



( -S12C23 S12S23 
C12C23 -C12S23 

\ S13 S23C13 C23C13 J 



M h M™ G = 

V 



C12C13 -C12S23S13 -C12C23S13 

S12C13 — Sl2 s 23 s 13 — S 12C23 S 13 





Mi. 



(30) 



For Model b(KM)), we have 



M, 



KM 
11 



C\ S1C2 S\S2 
S1C3 CiC 2 C 3 CiS 2 C 3 

V -S1S3 C1C2S3 C1S2S3 / 



M l , M, 



KM 
12 



(0 





) 







~S2S3 


C2S3 


Mi 


V° 


S2C3 


-C2C3 J 





(31) 



The Jarlskog parameters J and the effective mass < mpp > in the two models are the 
same as those in models a(PDG)) and b(KM)). 



III. NEUTRAL HIGGS AND CHARGED LEPTON YUKAWA COUPLINGS 

The different ways to identify the CP violating phase in the lepton sector with that from 
the spontaneous symmetry breaking sector also restrict the forms of the Yukawa couplings in 
the model differently. New interaction due to Higgs exchange can generate some interesting 
phenomena. Here we display the neutral Higgs coupling to the charged leptons for Models 
a) and b) and discuss some consequences. We have 
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Model a) : C llh o = -Z L M,Z fl 
Model b) : C Uh o = 



r vf 2 v s 



y (H° 2 -ia 2 ) + ±H°+ V -^(H°-ia) 
N A v 3 v N a 



+ h.c. , (32) 



rj- * 



Mi + ^VpMNsMn)l R {H^ 



%a\ l 



+l L M l l R 



I 

N a 



+ h.c, (33) 



N A - ~ v 

where I^l shown in above two equations are charged lepton mass eigenstates. Other Yukawa 
couplings are given in the Appendix. 

It is interesting to note that for Model a), the couplings are flavor conserving, but for 
Model b) the couplings can have neutral flavor changing current at the tree level. This can 
be used to distinguish these two models by looking at, for example \i —> eee and r — > I1I2I3 
decays. These decays have not been observed, but there are stringent bounds from various 
experimental measurements with 



Br(r -»■ eee) exp < 3.6 x KT 8 [24 
Br(r ->■ /uee) exp < 2.7 x 10 
Br(r ->■ /i/ie) exp < 2.3 x 10 



24 



24 



Br(r -»■ wfl) exp < 3.2 x 10 -8 
Br(r ->■ e/ip) exp < 3.7 x 10 -8 



24] 



25] 



Br(/x -)■ eee) exp < 1.0 x 10" 12 26]. 



Br(r -> ee//) exp < 2.0 x 10" 8 [24] 



(34) 



In the following we study these decay modes in more details to see if it is also possible 
to distinguish the Model b(PDG)) and Model b(KM)). We have 



PDG M PDG 



PMNS JKJ /1 



'13 



S23 s 13 c 13 



C23 s 13 c 13 



S23S13C13 C23 + S23C13 



-S23C23S13 



Mi. 



(35) 



V C23C13S13 -S2 3 C23S 2 3 sj 3 + cl 3 C 2 13 J 



1/KM M 
V PMNS IV1 11 



KM 



1 



C2 c 2 s 2 



Mi 



(36) 



\0 c 2 s 2 si J 

These two models have different predictions. For Model b(PDG)), at the tree level it 
is possible to have all the 7 decay modes listed in Eq. ( 134"|) . but for Model b(KM)) it only 
allows r — > fMfMp, and r — >■ fiee to happen at tree level. The decay rates of these processes 
depend on the PMNS mixing angles and phase, the Higgs masses, and also the two VEV's 
V\ and v 2 . To have some concrete feeling about how experimental data put constraints on 
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PDG parametrization KM parametrization 




02468 10 02468 10 

v, v, 

FIG. 2: The contour plots of v\ and t>2 in the range of 1 to 10 GeV for the upper bounds of fj, — > eee, 
and r — > hfoh- Here li can be one of the e or /i. The directions pointed by the arrows are regions 
satisfying the upper bounds. 

the models, we take the PMNS parameters with their central values described earlier, and 
the Higgs mass mediating FCNC interaction to be 100 GeV for illustrations. With these 
parameters fixed, the branching ratios depend on v\ and v%- Given the upper bound of a 
branching ratio, one can finds the allowed regions of 

We show the results in Fig. 2. We see, from the figure for Model b(PDG)), that the 
combined upper bounds fi — >■ eee and r — >■ e/i/2 give the most severe constraints requiring 
v\ to be larger than 7.4 GeV and t>2 to be larger than 2.5 GeV. Once these are satisfied, 
the other bounds are satisfied. For Model b(KM)), r — > gives the strongest constraint. 
The allowed region for for Model b(PDG) is bigger than that for Model b(KM)). For 
Model b(PDG)) it is possible to have /i —> eee to be close to its upper bound, but for Model 
b(KM)) it is zero. Future experimental data can distinguish these two types of models. 

IV. CONCLUSION 

We have studied some implications of models for identifying the spontaneous CP violating 
phase to the Dirac phase in PMNS matrix. This identification forces the Yukawa couplings 
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for lepton sector only dependent on lepton masses, mixing angles in PMNS matrix, and 
VEV's of Higgs bosons. Using different parameterizations of the PMNS matrix results 
in different models. We have studied in detail two popular, the PDG and the original 
KM, parameterizations. We then studied effects on properties of the PMNS matrix and 
neutrino masses. Since the CP violating phase in the CKM mixing matrix is determined by 
experimental data, CP violating phase in the lepton sector is also fixed. The mass matrix 
for neutrinos is constrained leading to constraints on the Jarlskog CP violating parameter J, 
and the effective mass < rripp > for neutrinoless double beta decay. The Yukawa couplings 
are also constrained. Different ways of identifying the phases have different predictions for 
ji —> eee and r — > lil 2 T 3 . Future experimental data can be used to distinguish different 
models. 

Acknowledgments This work was supported in part by the 985 project, NSC and NCTS. 



Appendix A 



To obtain the Yukawa couplings, it is best to work in the basis where un-physical Higgs 
have been removed. The un-physical Higgs bosons are the Goldstone fields h w and h z "eaten" 
by W and Z. There is also the axion field a which is invisible in this model. They are [lo| 

1 



h 7 



(vi^ + v 2 h 2 + v 3 h 3 ) , 
(vxAx + v 2 A 2 + v 3 A 3 ) , 



(-ViV^Ai - V 2 v\A 2 + v\ 2 V 3 A 3 - v 2 v s A s )/N a 



(Al) 



where v 2 = v f + v \ + v\ and A 7 " 2 = v 2 [y\ 2 v\ + v 2 v 2 ) with v 2 
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The physical fields, ai j2 , a and Hf related to the original fields are given by [10] 



A 2 
A 3 
K A sJ 

h 2 
\h 3 ) 



V2/vn 


-Viv 3 v s /N A 


vi/v 


-viv 2 /N a S 


' a\ 


-vi/v 12 


-v 2 v 3 v s /N A 


v 2 /v 


-v 2 v 2 /N a 


a 2 





vf 2 v s /N A 


v 3 /v 


vl 2 v 3 /N a 


K 





vf 2 v 3 /N A 





-v 2 v s /N a ) 


V a J 



v 2 /v 12 v x v 3 jvv l2 v x jv 



\ 



-vjv 12 v 2 v 3 /vv l2 v 2 /v 
-v 12 /v v 3 /v J 



(ht\ 

H 2 
\h w ) 



(A2) 
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where N\ 



| + v 2 v 2 ). d\ t 2 and H 12 are the physical degrees of freedom for the 



Higgs fields. With the same rotation as that for the neutral pseudoscalar, the neutral scalar 
Higgs fields {Ri, R 2 , R 3 , R S ) T become {H® , H 2 , H® , H®) T . Since the invisible axion scale v s 
is much larger than the electroweak scale, to a very good approximation, N a = v 2 v s and 
N A = v 12 vv s . 

Rewriting the Yukawa couplings in terms of the Higgs fields shown above, Hf, di and 
axion a, we have 



'—— — + —M m ) N R H{ 

V 2 V 12 V\V 2 



iV2 f l L —M D N R + l R — M x v L ) H 2 + h.c. 
vv 12 vv 3 



(A3) 



The couplings of neutral Higgs with neutrinos are given by 

£>vvh° = — W, 



-M D — + M m ^-) (#° + ia x ) - ^M D (# 2 ° + ia 2 ) 

V2V12 ViVn/ Na 



Md u0 v*M Dfv0 



V 



V2V12 V1V2 
{Hi + ia 



N R --{N R yM N N R 



v 



{H% + ia) 



+ h.c. 



(A4) 



For Model b), Yukawa couplings of charged Higgs and neutral Higgs with neutrinos are 
shown as 



^-vlh ± 



vv 3 



1- / V\ 

l VpmnsM d N r H 2 -\J2iv L 



Ml + ^-M n ) l R Hf 

V 2 V 12 v x v 2 



fn — V3 



vv 12 



M\l R H 2 + h.c 



C 



uuh° 



-VZM D N R 

1 



V ^{H» + ^a2) + -Hl + V -^{Hl + ^a) 

N A v 3 v N a 



R 



r ^3 m o ■ N 

[H 2 - ia 2 ) 



-{Hl-ia) 



(A5) 



(A6) 



-AW * " N a 

Note that in Eq. (1A3HA6j) II and l R are already as their mass eigenstates. The neutral Higgs 
couplings to charged leptons have been given in the text. 
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